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ABSTRACT 


On the evening of the 21 October 1977 an experiment 
was run in a cross-section of an east-west segment of the 
North Saskatchewan River Valley in Edmonton Alberta. Mea- 
surements of temperature, wind, and carbon monoxide concentra- 
tion were made in an attempt to determine the micrometeoro- 
logy of the valley. 

The inversion in the valley was more intense than in 
Eheseityawhichwasesi Cuatedyon the plained, The south tside of 
the valley was colder than the north side, and its tempera- 
tures were comparable to a rural station. 

Drainage winds were experienced on the slopes of the 
valley. The along-valley wind shifted from easterly to 
westerly as the inversion formed. During this time the 
prevailing wind remained southeasterly. A 15 minute periodi- 
city in the drainage wind on the south side of the valley 
was determined by spectral analysis. 

Carbon monoxide concentrations were found to be 
strongly influenced by the drainage wind, with a peak in 
concentration three hours after the late afternoon rush- 
hour traffic maximum. 

The data were not of sufficient resolution to describe 
the complete valley wind circulation, but the potential for 


pollutant entrapment within the valley was established. 
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CHAPTER 1] 


INTRODUCTION 


Wot) Reasons (for the? Study 

Many cities the world over are built around smal] 
river systems, where the river valley is below the level 
of the city. In many cases, these valleys are used for 
public recreation, including parks, amusement areas and 
recreational boating. These valleys are then at times the 
lOcavrOn ss Oyeslandemcongregaltons s0tspeOvle a WihlCh ein .vurn 
leads to high levels of motor vehicle activity in access 
roads and parking lots. Those areas of the valley not used 
for recreation are often sites for major urban roads and 
freeways, as well as various types of commerical enter- 
prise. Due to the small size of these valleys, the 
concentration of pollutants would not be expected to be 
of concern during the day when thermal mixing processes 
transport the pollutants out of the valley; but what of 
the night? Especially during the summer, evenings are 
when the maximum concentration of people would be expected 
and when the formation of inversion conditions would be 
most likely to concentrate the pollutants in the river 
valley. The questions arise: will the inversion conditions 
formed in these small river valleys be independent of the 


larger urban environment?; how will the urban environment 
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affect the river valley?; will the circulation within the 
valley be a closed system? There seems to be few studies 
that have concerned themselves with the micrometeorology 
of small urban valleys, but perhaps some of the observed 
circulations in large mountain and river valleys would 
apply to the small river valleys. 

The phenomenon of mountain and valley winds has been 
investigated by Geiger (1965), Defant (1951) and others. 
Observation shows that well-developed local circulations 
with marked diurnal variations are formed in valleys lead- 
ing into mountain ranges. There is a double system of 
periodic winds; during the day there are upslope and 
upvalley winds and at night downslope winds and downvalley 
winds. The slope winds are always initiated first and the 
valley winds follow. An idealized description of valley 
Winds iS given by ™Defant (1951); (Figs 1.1) At sunrise 
upslope winds begin (white arrows), but the valley is still 
colder than the plains and the downvalley winds (black 
arrows) are still blowing, fed by return circulation from 
the slopes. It soon dies out with further heating. In 
the forenoon (B), the slope winds are strong and there is 
a transition from a mountain wind to a valley wind, with 
the valley and plains temperatures the same. By noon and 
early afternoon (C) the slope winds are diminishing and 
there is a fully-developed valley ance with the valley 
warmer than the plains. In the late afternoon (D) the 


slope winds have ceased and the valley wind continues. In 
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the evening (E) downslope winds begin and the valley winds 
diminish. The valley at this time is only slightly warmer 
than the plains. By early night (F) there are well- 
developed downslope winds and a transition from a valley 
wind to a mountain wind. The valley and plains are at the 
same temperature. The middle of the night (G) sees the 
mountain wind fully developed, with the valley colder than 
the plains and by late night (H) the downslope wind has 
ceased and the mountain wind fills the valley. This 
description of the progression of the valley winds will 

be further complicated if there is a cross-valley prevail- 
ing wind.(Tang; 1975), which may result in closed cell (cut 
Of P)yweircula tions amuchis then iformsetheebas is on whitch tte 
begin the investigation of small urban river valley circula- 


tion and gain an insight into possible pollution problems. 


he2gsednomton iseRiver Vahey 
Within athe tity cof eEdmontometss 230 UNS iisl3OnWy the 


dominant topographical feature is the North Saskatchewan 
River Valley. A large part of this valley is used for 
public recreational activities and includes much of the 
city siipankland. gaboreth is areason aire isi tof winterest wto 
study the valley microclimate and pollution potential and 
if possible develop simulation (computer) models. Within 
the valley are complex arrangements of motor vehicle pollu- 
tion sources with bridges (cross-valley line sources), 


streets and roads (along-valley line sources) and parks 
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and residential districts (area sources). This valley is 
Known to possess its own microclimate, with, for example, 
minimum observed temperatures on clear nights on the floor 

of the valley comparable to those reported at nearby rural 
Stations (Hage, 1972). Evidence for poor mixing within 

the river valley is occasionally seen in fog situations 

where the layered nature of the fog is indicative of the 
absence of effective mixing between valley and plains air. 
Klassen (1962) found that the river valley system in 

Edmonton simulates the description of the model suggested 

by Defant (1951), but he did not observe an upvalley wind 
which was independent of the prevailing airflow. He observed 
a downslope wind under conditions of clear skies and light 
prevailing ince found the downslope flow to be weak, 

very unsteady and quite shallow. Air was seen to flow down 
the valley even when the prevailing wind blew perpendicularly 
across the valley. It was suggested that a double-vortex 
circulation develops under favorable conditions, with 

ascent over the relatively warm river water and descent 

along the slopes. 

In an attempt to answer some questions on the circula- 
tion within the valley, field experiments were run within 
the river valley in the center of Edmonton. Measurements 
of wind, temperature and carbon monoxide concentrations 


were made. 


Voc '? an? oo sidetn ne 


oe phase 
ae 


a? uti fas ehh! 
off Wak +o? Aste 


aayy “oinen Jy batseain oid ie ; 
ER a | fae 
sotta gn me rises ert | 

w hyagthen-ee ‘wer alee bavay 
beg ol tew: pe snuna a ae 
java “ol Pan Aqwant ‘ait tenn, brag) yaar, nese 
fohas ont Fo no higes ily ony sit efumte 0! 
facny civ ait dea soncSety ee ged’ , (120) Paeret 
’ Ctewhng GTS rnobwsqeha?, view ltd my 
: he’. and? (Sher Sebng. bat ¥ cant apn @ 
eo ft jy ate beve® bene bety ont ao 
noe@ atworth elt Oe WT: bas ybestanu vray : 
+e weta bai ont th ewwere add Hetlw nova’ wert ay od 
;' ful , duode bot eteeeelegw th . Vat Tavcees ragtee 
( Seo TS TRACI eee ei agai evee maga 


suganah bap Wedsae 16+) ee wphiytanten stt “540° $nases. 


oe: “ wanste ant booth oe 


Tale i, wei? Ang EVohe eanh mag +~nMeie 6) aqme dm as<atl i.  @ 


iphw voy eae dimebienee otal ayet ley ens maid - 
J 

ast uT ee “e68n08bE) te ‘knee od? wi vat Tay Shay 

soi ¢47 teeaneD Sebi scned nodena tae everett : ve 


i 
7 a 


- 


me 


CHAPTER 2 


EXPERIMENTAL LOCATION AND INSTRUMENTATION 


2.l4 Lopogrdapny 
The City of Edmonton is located on the banks of the 


North Saskatchewan River in central Alberta. Edmonton lies 
670 m above sea level on a plain that slopes gently from 
the southwest down to the northeast (Fig. 2.1). The North 
Saskatchewan River meanders through Edmonton from the south- 
west to the northeast and the valley is steep sided and 
narrow, with an average depth of 50 m. The width of the 
valley varies from 1000 m in the straighter sections to 
1500 m. Data were collected from the Edmonton City Tower 
(A), the Municipal Airport (YXD), the International Airport 
(YEG) and mainly from the High Level Bridge (B) area. The 
International Airport is situated 32 km due south of 
Edmonton in a rural setting. 

The detailed study of the river valley was done in 
a cross section of the valley in the vicinity of the High 
Level Bridge (Fig. 2.2). The bridge consists of two traffic 
levels with a railbed on top and a motor vehicle roadway 
below. Pedestrian walkways are found on both sides of the 
roadway. The walkways are 42 m above the river. The south 
slope of the river valley in the area of the bridge has no 


major roads running parallel to the river and no through 
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roads west of the bridge below the rim of the valley. The 
south slope is heavily treed. The north slope vegetation 
consists of tall grasses and is sloped more steeply than 
the south slope. River Road is a major traffic artery 
within the city. The road down the north slope passing 
under the, bridge as lightly travelled) Trataric counts 
along River Road vary from 1,700 motor vehicles per hour 
ate? /O008MDT tol 200 per hourvat O100,;MOR. Lhe south rim is 
level and there are large high-rise buildings running for 
most of Its length in the vicinity of ‘the High Level Bridge. 
The north rim of the valley is a large open area of short 
grass, with no nearby large buildings. 

The locations of the experimental stations are indi- 
Catedi@ins fig. 2.2, The numbers 1 to 5 indicate the loca- 
tions of the main fixed experimental stations and the 
letters A to H indicate the locations of the carbon monoxide 
monitoring stations. Stations 64and © are noteindicated 
as both were located beyond the area covered in Fig. 2.2. 
Station 6 was located on the south rim of the valley about 
0.5 km west of the bridge and station C was located at the 
base of the north slope WNW of station D. The elevations 
of the main stations (1-6) are indicated in the cross sec- 
tion of the river valley 40 m west of the High Level Bridge 
(CFageics3)s so tateons: 4a-to 4d) aindicatesithe locations of 
a thermocouple array that was suspended from the High Level 


Bridge. 
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During the summer and fall of 1977, experiments were 
run to provide data for modelling studies. Temperature, 
wind and carbon monoxide concentrations were measured under 
light winds and clear skies. Thus the microclimate of the 
valley was observed under inversion conditions. The first 
experiment, run on July 13, was a preliminary study to gain 
familiarity with equipment and to establish procedures. 

The second experiment, run on July 20, did not yield light 
winds, and the mean wind penetrated to the floor of the 
valley. The third experiment held on July 26 gave the 
required conditions. 

It became obvious that more was occurring in the 
micrometeorology of this small valley than had been expected 
and, in many respects, the findings resembled those for 
larger valleys. A survey of the literature did not reveal 
any work done on this scale of river valley, although many 
studies had concerned themselves with mountain valleys or 
large river valley systems. 

To investigate the observed complexity, additional 
instrumentation was obtained or designed to allow study at 
locations where only battery power could be used. As very 
light winds had been observed previously, this equipment 
also had to have good sensitivity and stability. 

A fourth experiment was run an October 21, during the 
period of 1800-2400 MDT, although some stations have slightly 


longer records. Additional information was gained from the 
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Municipal Airport, the International Airport and the Edmonton 


Gity Tower. 


2.3 Station Instrumentation for Experiment 4 


The positions available for the stations were all 
easily accessible by the public, as much of the river valley 
is part of the regional parks system. As a result all sta- 
tions were manned continuously and all equipment was removed 
at the end of each experiment. Portable meters - were used 
for the measurement of carbon monoxide concentrations and 
the output recorded on portable chart recorders.* The meters 
and recorders were carried from station to station and five 
minute readings of concentration were taken at each station. 
The instrumentation which was located at each station is 
listed in Table 2.1. All thermographs were held in Stevenson 
Screens, with those at stations 1 and 3 suspended above the 
ground, those at stations 2 and 6 placed on the ground and 
the screen at station 5 was placed on a picnic table. Tempe- 
rature profiles were obtained at station 4 by suspending 
thermocouples from the bridge at 5, 13, 21 and 29 m below 
the road level of the bridge. The anemometers at stations 
1 and 3 were supported by standard television antenna tower 


sections and the profiling anemometers at station 2 were 


Interscan Corporation model no. 1148. 


Esterline Angus Corporation model no. MSH401BB. 
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Table 2.1. Meteorological Instrumentation 


Station Variable Instruments 

] wind speed and direction cup and vane anemometer 
air temperature thermograph 

2 slope wind propeller anemometer 
wind profile cup anemometer 
air temperature thermograph 

3 3 wind components propeller anemometers 
wind profile cup anemometers 
air temperature thermograph 
water temperature thermograph 

4 bridge level air temperature thermister 
4 level air temperature thermocouples 

5 slope wind speed propeller anemometer 
wind direction vane 
air temperature thermograph 


6 air temperature thermograph 


a = ' ee eee 
ithat f SisyY. O48 | | 
it 
r rn anor CC . “bits bee 


Se uinese aa were athe’ in ua 


lnevvorea? a prude rage ate = 


: i |. fi 
ane ash fegowe: TeeesOGmO? bn ¢ 
i¥ = AN Banh 
us ay % ef Movg bate ae ny 
‘rd sr aqmes its” te 
wit? asu7a19omel Tedaw 7 
| g ee a 


sod 9 )erertd oy Sede ata (see! oles ee! 


‘ony arts aruserveqse? Vial overt». 7 oseip = 
vk us ears 


rome ne 15) hones hwets bitte aqole an , 


-; : i _ ig hs b 
one ae ies anh, aay 


Aqsa Tyce pete tabi) 7 Sy 3 


igeteomiges 


9 


ee ne © 


14 


Supported by a guyed metal tube. The propeller anemometers 
at the river valley slope stations 2 and 5 were aligned along 
the slope so as to measure the up and down slope component 

of the wind. The heights above ground of all instruments 


ane listed in Table 2:2. 


2.4 Calibration 

The thermographs used were made by Casella, Thies and 
Ryan Instruments, | the latter being submersible for obtain- 
ing water temperatures. Measurements of temperature at sta- 
tion 4 were obtained from a thermister thermometer. Tas 
thermometer was found to be accurate when calibrated against 
a standard mercury-in-glass thermometer and was used as a 
reference for all other instruments because of its rugged 
portability. Profile measurements were made by thermocouples 
of nickel-constantan of 0.013 mm diameter. This size was 
used to avoid radiation effects, as the elements could not 
be shielded. The output from the thermocouples was measured 
by a digital readout> which was not capable of reading below 
a machine reading of O°C. Calibration of the temperature 
measuring instruments was carried out at several tempera- 


tures, using both the Meteorological Division cold room and 


Ryan Instruments model no. F-8. 


2 Weather Measure Corporation model no. T618-1-B. 


Omega Engineering model no. 200-EC2-DSS. 
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Table 2.2. Height of instruments above ground. 


Instrument 


Station Height (m) 


thermograph 


cup anemometers 


propeller anemometers 


Vanes 
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still outside air, over a temperature range of +20°C to 
-15°C. Many of the instruments used had slight temperature 
dependencies and for these regression lines were fitted to 
the calibration points. Accuracy of the thermographs was 
affected by two factors: the extraction of data from the 
charts, and the slightly differing positions of the paper 
on the recorder drum. Data could be extracted from the 
chart itself to within 0.1°C and the positioning error can 
be inferred from the average standard error of regression, 
which was about 0.2°C, implying a position error of about 
0.1°C. The thermister thermometer could be read to within 
0.1°C and the digital readout for the thermocouples is 
specified to within 0.1°C resolution. 

Impulse counting anemometers were used to obtain wind 
profiles and a cup anemometer and chart recorder system? was 
used at station 1 to measure prevailing wind. All pro- 
pellers were of the Gill type. Anemometer calibration was 
carried out in the University of Alberta Mechanical Engineer- 
ing Department low speed wind tunnel, which is capable of 


ie The calibration reference for 


] 


speeds down to 0.04 ms 
the speed range of 0.05 to 2.0 ms was a hot film anemo- 
meter owned by the Mechanical Engineering Department and 


which was calibrated from first principles and linearized 


Atte, MSer en so Model AST. 


° R.M. White anemometer model no. 5-902, recorder 
model no. 1-188. 
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1 


over this range. Above 2 ms the Gill propeller anemo- 


meters were used as a reference as the response was found 
] 


to be linear above 1 ms The’stallespeed of thesRemM. 


White cup was found to be 1.5 m ae The Gill anemometers 


were very close to one another in calibration and stall 
] 


speed, the latter varying from 0.15 to 0.18 ms The 


counting cup anemometers also proved reasonably sensitive 
with stall speeds grouped about 0.32 m cant 

The zero calibration of the carbon monoxide meters was 
found to be unreliable, so to detect zero drift, the meters 
were left to run continuously for the duration of the experi- 


ments «In this way the zero drift could be observed and 


corrected for, as the drift was relatively slow and steady. 


ous GChart Recorderabdesign 


During analysis of Experiment 3, it became obvious 
that a sensitive chart recorder was needed to obtain satis- 
factory resolution of wind speed, which could change by an 
order of magnitude or more during the period of interest. 
Due to site location, the recorders had to be light and 
battery powered, since none of the sites had access to 
commercial power lines. As Rustrak® chart recorders were 
available, it was decided to design electronics for these 
recorders. The Rustrak recorders consist of a one milli- 


ampere meter movement, whose indicating arm can be clamped 


6 Rustrak Instrument Company model no. 214C. 
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by a striker bar to mark pressure sensitive chart paper. 

The chant spapemisedriven bya 12%volted.c. motor, which 
alsomcontrors mine, striker bar. 2 lhe tdrdin of athis motor); 1s 
very low, so that the current drain on the 12 volt lantern 
battery used will not reduce the output voltage by a measur- 
able amount by the end of the experiment. Depending on the 


] 


gear train selected, speeds of up to 30.5 cm hr (12 inches 


hr7!y could befobtained. “During Experiment 4,-station 5 


] 


used sa chart speed of 15.25 cm hre — (6) inches Are) and all 


other stations used 30.5 cm hr’, 

Operational amplifiers were selected for the chart 
recorder electronics due to their precisely-defined para- 
meters and ease of replacement in the event of failure. 
Originally only 741 type operational amplifiers were used, 
but due to their high input bias requirements of 2 x 1077 
amperes, any change in the input or feedback loop impedance 
caused large changes in the output offset. No rearrangement 
could cure this problem due to the change of three orders 
of magnitude in input impedance shown by the required input 
attenuater. After some experimentation it was decided to 
use the National Semiconducter LM310N as an input buffer. 
hi Stamp limi er ghias “aeitypical: input bias current of.) x 107? 
amperes and isolates the 741 from impedance changes. The 
current feedback loop around the 741 (Fig. 2.4) was found 
to provide the most stable operating point, with no instabi- 


lities developing and with no discernible output zero drift 


after several hours. 
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As shown in the circuit diagram, SW] is the input 
attenuater, stepped in a semi-logarathmic sequence and pro- 
viding a constant input impedance of 50 K®, as the shunt 


impedance of the LM310N is 1.5 x it 


ohmsas EAdleyesistors 
in the attenuater are one percent metal film resistors for 
best stability and accurate input attenuation. The large 
range of input voltages available, while unnecessary for 
Experiment 4, was designed to allow greater flexibility in 
future experiments. Resistor R3 limits the maximum current 
available to the chart recorder, should a device failure 
occur, to 4 milliamperes to prevent burning out the meter 
movement of the recorder. Output offset is controlled by 
Variable resistor R2s in conjunction with the internal cir- 
CUILTEY OT (CNC 7412 Gain 15 set bY resastors) R4 and) RS, with 
the maximum gain limited by R4 to two and the minimum gain 
limited by R4 and R5 in series to 0.7. Switch SW2 is the 
power on-off switch. There is interaction between the out- 
put offset control and the gain control. Changing the gain 
will change the offset, but the offset will not affect the 
gain. The initial set-up procedure consists of adjusting 
both gain and offset one after the other until the full- 
scale range is accurate. Once the full-scale reading is 
correct, the offset may be adjusted and the gain will not 
be affected. Gain stability depends on the stability of 
resistors R4 and R5 and output offset Stability will depend 
Onstne stabaiity of Re. Ihe cincuit is: not voltage sensi- 


tive and will operate well at 7.5 volts. Total current draw 
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Of theca rcunteat, tullesicalverde fection its 5.3 milliamperes 
from the positive supply rail and 4.3 milliamperes from the 
negative supply rail. 

The performance of this design was tested on all 
ranges. From maximum to minimum gain, the zero offset 
Mart essby, OL Se percent: and-fuliimscale accuracy. isewitthin 1.0 
percent. Total error of any reading was 2 percent with all 
errors additive. If a limited number of closely spaced 
ranges were used, accuracy was typically within the 1 per- 
cent level. The circuit was also tested with a calibrated 
voltage source for temperature drift. Within the range of 
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CHAPTER 3 


RESULTS OF TEMPERATURE MEASUREMENTS 


3.1 Temperature Profiles 


The temperature profiles obtained at station 4 were 
augmented by temperatures measured at stations 1] and 3. 
The assumption was made that at any given height in the 
valley, excluding the slopes, the temperatures would be 
uniform over the horizontal distances separating the sta- 
tions. To extend the temperature profiles into the urban 
environment, temperature measurements from the 10 m (A), 
Ah? Vic B anc! ole oem Gc) levels one the Edmonton City lower 
were used (Fig. 2.1). The actual temperatures could not 
be used as comparison of the tower temperatures with 
temperatures obtained from the municipal airport (YXD) 
during inversion and lapse conditions showed differences 
of up to two degrees between the two locations. That these 
temperature differences were due to location was demon-~ 
strated when periods of continuous rain were used for 
comparison. During such occurrences, temperatures over 
large areas can be expected to be uniform and in these 
cases the tower and YXD temperatures were in good agreement. 
The assumption was made that the shape of the temperature 
profile would be the same over the river valley as at the 


location of the Edmonton City Tower and that the shape of 
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the temperature profile would be maintained with respect 
to its height above ground. To relate the tower tempera- 
tures to temperatures measured at the river valley, an 
estimate of screen-level temperature expected from the 
tower profile was made by extrapolating the 10 m to 45.7 m 
temperature difference to screen-level. The expected screen- 
level temperature was then equated to the actual tempera- 
ture at station 1 and the entire tower temperature profile 
adjusted by the amount of the difference obtained. The 
tower measurements of temperature were assigned to the same 
height above ground at station 1 as at the Edmonton City 
Tower location. Temperatures obtained at half hour inter- 
vals, including the altered tower and river water tempera- 
tures a@reailisted in Table 3.1. 

Temperature profiles plotted every hour from 1800 to 
2400 MDT are shown in Fig. 3.1 and, as can be seen, the 
variation of the profiles with time is reasonable. At 1800 
the profile within the valley was isothermal and temperatures 
decreased with height above the valley. This is reasonable 
as the valley at this time had been in shade longer than 
the air above the city. A peak in the temperature profiles 
at bridge level in this and other profiles was attributed 
to the thermal effect of the High Level Bridge and the 
traffic on it, as the thermometer was located on the down- 
wind side of the bridge. The bridge level temperature was 
much larger than expected from the other profile tempera- 


tures (up to 1.5°C), as well as erratic, so it was not 
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included in the profiles. The thermometer was moved after 
2100 to try and minimize the thermal effects, but the 
temperatures were still unreliable. 

At 1900 cooling was evident in the lower 25 m of the 
river valley and the temperature profile above the city was 
isothermal. By 2000 the cooling in the lower levels of the 
valley was very pronounced and in later profiles this cool- 
ing had extended upward. Once the inversion had formed, 
the temperatures at the bottom of the river valley were 
always lower than would be expected from the Edmonton City 
Tower profile, indicating that the river valley may have 


been closer to a rural than to an urban temperature regime. 


o.ce TIMeaniablons 
In order to follow the variation of temperature with 
time, the temperatures were plotted at ten minute intervals 
with the axis of one station plot shifted from another to 
illustrate the simultaneous time variations of temperature 
at all stations. The station temperatures used for the 
profile measurements within the valley are shown in Fig. 3.2. 
As can be seen, the onset of the inversion at station 3 was 
abrupt, starting after 1940 MDT. The temperature fell 
rapidly to a plateau at 2100. Rapid falls followed by a 
leveling off of the temperatures occurred three times before 
2400. The inversion can be seen propagating upward through 
Stanrens 4d tomy. mareed DY Che stick Marks In bid. soc, at 
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ture between 2100 and 2110 at station 4 was the result of 
moving the thermometer in an attempt to avoid the thermal 
effect of the High Level Bridge and was not an actual 
temperature decrease. The most probable real temperature 
trend is indicated by the dotted line. The initial 
accelerated drop in temperature at station 3 and the sub- 
sequent plateau was reflected in the temperatures of all 
other stations excep& station 1, tndteating that this 
phenomenon was a characteristic of the river valley. The 
mechanism responsible for this behavior is at present 
unknown. The second rapid fall in temperature after the 
first plateau was again reflected in levels 4d to 4a, 
aithough not as clearly asthe first décrease. The third 
rapid temperature fall at 2120 at station 3 was of greater 
magnitude than the second and was reflected in a rapid fall 
in temperature at station 4 one hour later than at station 
3. This is in agreement with the upward propagation rate 
of the initial inversion which occurred at 1940. 

The; two slope stations 2 and 5, the south rim sta- 
tion 6, the municipal airport (YXD) and the international 
airport (YEG) temperatures are shown in Fig. 3.3 in the 
same tormpeas, in Fig. 3.2. Statdons 5 and 6 showed no 
periods of rapid temperature decrease, but exhibited a 
steady cooling rate. The airport temperature showed large 
temperature changes from hour to hour, but the linear trend 
obtained with a least squares fit and indicated by a dotted 


line, showed general agreement with the cooling rates of 
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the other stations. Station 2 did show one period of rapid 
COOMIMGeDeDuinInGg at 2010. mndicated  pyrattick mark... Ihis 
time is in agreement with the beginning of the onset of 

the initial inversion shown at station 4c which was close 

to the height of station 2. The absence of further rapid 
decreases in temperature is probably due to the formation 

of a downslope wind, isolating station 2 from direct inter- 
action with the valley air. The absence of any such initial 
temperature drop at station 5 is surprising until the actual 
temperatures are studied (Table 3.2). The temperatures at 
station 5 were lower than those at any other river valley 
station and were comparable with those at the rural station 
YEG. The air temperatures above the rim of the valley at 
Station 1] were comparable to those at the urban station YXD 
as expected. The north slope and river station (2 and 3) 
temperatures were midway between urban and rural values. 

An explanation for the lower temperatures recorded at 
station 5 may be discerned by studying the change of tempera- 
ture with time from Experiment | (Fig. 3.4). As can be 
seen, the temperatures at both stations were comparable at 
1645, but station 5 began to cool earlier than station 1. 
When station 1 began to cool at about 2000, the rate of 
cooling at station 5 was similar and the temperatures at 
station 5 remained lower than those at station 1. The early 
cooling of station 5 is consistent with the early shading 

of the north-facing slope from the sun. This phenomenon 


was also observed during Experiments 2 and 3 in varying 
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Table 3.2. Comparison of temperatures (°C) at main river 
valley stations with the urban airport (YXD) 
and the rural airport (YEG). 


Station 

Time ] 2 3 b 6 YXD WEG 
1830 Cary 8.4 8.2 -- -- 

1900 [Ege 7.4 7.4 4.2 Sig 6.6 ‘beers 
1930 ewe Gag Gao 350 1B 38. 

2000 629 eno Sri0 Cig Ora Sad -.4 
2030 Gaz D3 aes 4 1.4 4.5 

2100 B22 3.4 209 05 4.4 es, Ba -.6 
2130 4.7 250 19 - Sy) 

2200 aiar EG 1.6 -.9 Siac 35:0 -3.4 
2230 SZ 1.4 Leet -1.4 ea 

2300 Saul pat -.4 -1.9 eel ifs -2.1 
2330 25 -.3 -.6 -2.5 0.6 
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degrees. Since this difference in temperatures appears to 
be common on clear evenings it suggests that an east-west 
segment of the valley is asymmetric in temperature and that 
little interchange of air occurs between the two sides of 
the valley, at least near the surface. This casts doubt 

on the existence of any type of double-vortex closed circu- 


lation within such segments of the valley. 
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CHAPTER 4 


RESULTS OF WIND MEASUREMENTS 


4.1 Wind Above the River Valley 

Wind speed at station 1, representing the above-valley 
flow, was measured by a three-cup anemometer and the out- 
put was recorded on a Rustrak chart recorder. As this 
chart recorder marks pressure sensitive paper with a 
StGMIKerEbDarnethe recorndingecons is tS Of a senves, Of dots. 
The density of the dots is too high to enable resolution 
of the time sequence of successive individual dO.US3 
although the trend was obvious. To extract the data, an 
averaging procedure was adopted, wherein the chart record- 
ing was divided into sections and the average value of the 
trace abstracted by visual observation. This procedure 
was tested for accuracy by abstracting values by the visual 
method, then averaging the same data where the value of 
every point was A erer ign adie Errors were found to be less 
LhaneOmpercent =i neal caseseand gless ethane o pencent tor 
most cases. The wind speed recording of station 1 was 


divided into two-minute periods by this method. 


The stall speed of the cup anomometer used was 1.5 ms 


For a wind speed of less than 1.5 m al. the indicated 


wind speed will be zero. Where the wind speed is below 


the stall speed for only part of the time, the average 
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value abstracted from the chart will be a measure of the 
percentage of time the wind speed was above 1.5 m ae 
Where the wind speed remained above the stall speed, the 
indicated average speed will be correct. 

As expected, the wind speed decreased with time and 
became light as the inversion formed by 2100 (Fig. 4.1) 
av Station 1.” fhe wind jdirection was from the SE for most 
of the time, except for the period from 2030 to 2230, when 
the wind direction became unsteady, backed to WNW for 35 
minutes, then veered back to SE. At the same time the wind 
speed dropped to low values while the wind was unsteady and 
picked up once the wind had established itself from a given 
direction. To see if this behavior was truly that of the 
mean wind, the wind directions from YXD were plotted, and 
the agreement was generally good considering the distance 
between the stations and considering that the winds at YXD 
are measured at the standard 10 m height above ground. 
Except for the directional shift to WNW at 2200, the mean 
difference in direction was 20 degrees with a standard 
deviation of 12 degrees. At 2200 the difference was 90 
degrees, indicating that this change in direction must 
have been related to the river valley and not to the pre- 
vailing flow. If the inversion, initiated in the river 
valley, is assumed to reach above the level of the rim of 
the valley, the time of arrival would be about 2050, which 
is the time at which the wind at station 1 began to back to 


WNW, indicating that the increase in wind speed and the 
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Prevailing wind direction and speed measured 
atestationed. «Wind)direction for municipal 
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Change in direction experienced at station 1 were related 
to the inversion propagating to this level. The tempera- 
ture at station, |. did. show anvincrease, in. the, rates of 
temperature decrease, but this was not any larger than the 
rate of temperature decrease observed at other times and 


Gannhote pe used to, contirm thesaryrival_of. the. anversion. 


4.2 Winds Within the River Valley 


nn ee 


facet ocathon. 3 

Wind measurements at station 3 were taken with a 
three-axis Gill propeller anemometer and two cup anemo- 
meters for wind profile measurements. Data were abstracted 
from chart recordings of the propeller anemometers by the 
Same method that was used for station 1, with averaging 
cymes of 170 to Taz minutes. the horizontal ly-positioned 
propellers were aligned along the valley axis (W to E) and 
across the valley (N to S). The wind speed and direction 
were obtained by assuming a cosine response for the pro- 
pellers and summing the individual results vectorially. 
The actual response of the Gill anemometers has been shown 
not, to pe a cosine (Drinkrow. 1972). but the error in this 
assumption is about 1] percent and is considered small when 
compared to other sources of error. The chart speeds were 
not constant over the length of the experiment, so to 
equate the two components, the chart speed was assumed to 


be constant between hourly time marks, and values of speed 
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were interpolated to the same time intervals of 1.2 minutes. 
The time series was then low-pass filtered to remove the 
high frequency components of the wind (periods of less than 
2.4 minutes), since these components are questionable due 

to the data extraction method and would not aid in the 
analysis. The trend of the time series was obtained by 
low-pass filtering over a one-hour time interval using a 
Bar he tact tmlbteanl.S.CemS.6.C,t1,0NimAn 2. uONeciS.cuss1on). 

Wind profile measurements were made at 1.4 m and 2.6 mm 
levels, with counting cup anemometers over periods of one 
hour. The average wind speed measured by these anemometers 
was found not to agree with the one-hour averages obtained 
from the propeller anemometers, in that the cup speeds were 
Significantly higher for low wind speeds. This effect is 
probably due to the wind speed being too close to the stall 
speed of the cups and to the resultant non-linearities. 

This means that, while the shape of the profile obtained 

with the cups will be good, the measured wind speeds cannot 
be equated to other sensors if the wind speed drops below 

the cup stall speed for a portion of the time. When the 
average wind speeds were greater than 1.0 m aa Cea Otel eS. 
including the propeller anemometer data, were reasonable 
(Fig. 4.2). After 1900, when the average speeds were below 
Tom Sa), ENCaswiCluS vO. O1rm=ehemp Gope.l ber data 1s dubious. 

The cup anemometers do, however, indicate that the wind 

speed increased with height before 2000 and became almost 


constant with height thereafter. 
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The wind speed at station 3 dropped to low values at 
2000, “Chen prekedsup= cor a constant speed of about,0.6 m 57! 
anter 2sO00n (hide 6.3).8 This 1s surprising, asethe wind 
speed would be expected to become very low with a well- 
developed inversion in the valley, but an examination of 
wind direction data indicates the reason for this (Fig. 
4.4). Before the formation of the inversion after 1940, 
the wind was blowing from the east and it began to veer at 
1940 to become westerly by 2040. The period 1940 to 2040 
was the time required for the inversion to completely fill 
the river valley. Thereafter the wind was always westerly. 
The westerly wind in the valley within the inversion layer 
was in accord with the change of wind direction to westerly 
Gib set debt Ola, lee This suggests the idea that the inversion 
propagated above the rim of the valley. 

There are several possible explanations for the 
Switch to a westerly wind. The river runs at a consider- 
able velocity in an easterly direction and frictional drag 
might pull the air with it. However westerly winds in the 
river valley with a southerly prevailing wind have been 
observed with the river frozen over. The slope of the 
river valley may also cause a drainage effect, but the 
Slope of the river is only 0.9 m/km and this does not seem 
to be a reasonable explanation for the rapid shift in 
direction in synchronization with the formation of the 
inversion. The most plausible theory is tied to the geo- 


graphical alignment of the North Saskatchewan River Valley 
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system. The river is generally aligned SW to NE and con- 
sists of a series of meanders that are generally aligned 
Neto’ SiandeWe tor Ee (Fige(2.1)eerThis®can* be represented 


schema ticald ye (Figt14: 5) .waAlthoughi once! ther inversionsis 


Compensating 
Flow 


Flow forced by 
prevailing wind Experimental 
Site 


II 


Prevailing Wind 


Fig. 4.5. Schematic of North Saskatchewan River Valley 
system showing possible cause of westerly 
winds at experimental site. 


established, the valley air is expected to be isolated from 
the prevailing wind, if the prevailing wind blows parallel 
to the valley axis, the air in the valley will tend to move 
with the wind and there will be an increase in mass in the 
leeward end of the valley and a decrease in mass in the 


windward end. Since the river valley is arranged stair- 
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case fashion, a southerly wind will move the air in the 
N-S aligned segments of the valley northward, and the 
resulting pressure gradient across the E-W segments of the 
valley would establish a westerly flow within those seg- 
ments. This could happen rapidly and could account for 
the rapid directional change in the valley wind. This 
theory, however, appears to be only part of the explanation, 
for a westerly wind in the valley in winter when the river 
is frozen has been observed with calm prevailing winds as 
evidenced by vertically-rising plumes from sources within 
the city. 

Vertical velocity was also measured at station 3 and 
recorded on a Rustrak recorder. As mentioned previously, 
it was not possible to distinguish the time sequence of 
successive dots and, in the case of vertical velocity, 
when the signal switches rapidly from positive (upflow) to 
negative (downflow), the previously-used averaging proce- 
dure will not give any significant information about 
vertical velocity. In order to derive statistical informa- 
tion, the values of each recorder trace were abstracted 
fon Jaaguven ecime period.) In this casewa period. of five 
minutes was chosen. Four such periods were abstracted and 
the means and variances were calculated. The variance of 
the vertical velocity with time 1s shown in Fig. 4.6, the 
final reading being an estimate, as the signal levels were 
too low to enable the value of individual dots to be read. 


After this time the signal was below the noise level 
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of the recorder, indicating no movement of the anemometer. 
Measurable vertical velocity activity did not end at 1900, 
but the activity cannot be described in a statistical way, 
aseratereactiviuy coccurned “in tsivontebursits:, vast stress 

built up, then released. The first measurable activity 
occurred at 1940, with a brief upward burst of 0.4 m oe, 
Four more bursts of upward motion occurred in the next ten 
minutes, then at 2000 activity occurred for five minutes 


Uo After this: there was- no 


] 


with peak velocities of 0.9 ms. 
activity Until cldasownen a veryepriet burst or O-2 mrs. 
occurred for 20 seconds. Vertical velocities were observed 


again from 2200 to 2208 with peak velocities of 0.3 m 57! 


and at 2215 with peak velocities of 0.3 m s7). After this 
time no activity was recorded. 

As a check on the consistency of the data, stability 
was determined from the available information at station 3 
and from the lowest level of the temperature profile 4d. 
The stability parameter is Z/L where L is the Monin-Obukhov 
length. As L could not be determined due to the unavailabi- 
lity of surface heat flux measurements, the assumption was 
made that Z/L was equal to the gradient Richardson number, 
R. (Panofsky, 1973). This assumption may be questionable 
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where g = acceleration of gravity 


fap) 
i} 


potential temperature 
T = mean temperature 
and where the partial derivatives were evaluated for a 


logarithmic profile using (Paulson, 1970) 


where F is, for example, 6 or U and the derivative is deter- 


Vie The wind speeds were taken 


mined for the height (ZZ) 
from the cup anemometers at 1.4 and 2.6 m. Although the 
values measured by the wind profile system are questionable, 
the change in velocity with height will be valid, as the 
errors in the measurements are common and should cancel. 

The temperatures were evaluated from the thermograph tempera- 
ture at station 3 and the thermocouple temperature from 4d. 
The change of temperature with height between these two 
levels was assumed constant. As the stability is sensitive 
to temperature, the temperatures were averaged one an hour 

to eliminate random errors and to match the one-hour wind 
speed averages from the anemometers. Fig. 4.6 shows that, 


CONSTSCeENtewrthecne time» oreonset of the inversion’ at? 1940 


MDT, the stability increased rapidly. 
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Ame. PO tatVons, 2 

The upslope (north) and downslope (south) component 
of the wind at station 2 was measured with a propeller 
anemometer, but the cross-slope component was not measured 
due to equipment problems. The data were abstracted and 
filtered and the trend calculated by the same process as 
that used at station 3. The averaging period was 1.2 
minutes and the trend was calculated over one hour. As an 
earlier study by Klassen (1962) indicated that the drainage 
flow was weak and shallow, cup anemometers were arranged 
toOmeasurel tne: Winds Drotilesate lo 2) andes mee The: cup 
anemometers used were of the same type as those used at 
Stations J -andathe earlier discussion of their damitations 
and errors pertains as well to the profile measurements at 
Sait Ont 2c 

In the early period of the evening the wind was blow- 
ing upslope, either in response to the mean flow or indicat- 
ing that a daytime type of valley circulation was still 
oecurning -(Figne 4.7). Ase the eshect of nocturnal cooking 
made itself felt, the wind at station 2 decreased in 
intensity and at 1940 began blowing downslope, indicating 
the establishment of a drainage flow. At 2000 the wind 
again blew upslope continuing until 2030, when it once 
again began blowing downslope and established a drainage 
flow, exceptefion a few brief; periods from 2230 to 2300, 
when the temperature remained constant. Although no wind 


direction measurements were available, the flow was con- 
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firmed 2s downslope by the station operator for the latter 
part of the experiment. From the temperature record, the 
inversion was seen to arrive at station 2 about 2010, 
although the maximum rate of cooling was not experienced 
until 2030, when the wind at station 2 became downslope. 
The timing of the arrival of the inversion and the activity 
of the wind, suggest that the upslope wind that occurred 
after the initial downslope regime was established and the 
inversion were connected. The most likely explanation for 
this 1s that the inversion Wayer acts as a ’fleor’* to the 
penetration of the prevailing wind so, as the inversion 
layer thickens, the mean wind flow will be constrained to 

a shallower depth in the valley and Ste affect the slope 
winds. Without a measure of the cross-slope wind component, 
however, it is not possible to determine if the flow is 

SE (the prevailing flow) or W (the channeling effect which 
resulted in the westerly winds at station 3). 

The wind profiles reflect the normal increase in wind 
speed with height, until the period 2030-2100 MDT, and are 
close to the theoretical logarithmic shape expected in the 
constant-flux layer (Fig. 4.8). From the period 2030-2100 
onwards, after the drainage wind had established itself, 
the profiles indicate a constant wind speed with height, 
Showing that the drainage wind, although light, was much 
deeper than expected. Observation of a steam exhaust plume 
from the Glenora Club did not show evidence of the down- 


slope flow at 2100, but by 2350 the plume was observed to 
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drifting southward, indicating that the drainage flow had 
extended in depth, on the river flood plain, to the top 

Of the building at oem. Evidence that ther drainage flow 
reached the height of the propellor anemometer array at 
station 3 can be seen from the cross-valley wind component, 


which was consistently from the north after 2230. 


eZ oe Sb CHONIsS 

The upslope (south) and downslope (north) component 
of the wind was measured with a propeller anemometer and 
the wind direction with a vane. The data were abstracted 
and filtered, and the trend was calculated with the same 
parameters that were used for station 2. A plot of wind 
direction was not produced as the wind was downslope from 
1830, when measurement began, to 0100, the end of the 
experiment (Fig. 4.9), and the measured wind direction was 
180 degrees: +815 degrees Aans noted in Chapter 3, ‘the 
temperatures measured at station 5 were lower than those 
at any other station because cooling began first at this 
station. It was also noticed that the temperatures implied 
a-difrerent regime at station 5,5 and this is reinforced. by 
the wind measurements. The downslope wind with no signifi- 
cant time variation in the trend implies a well-developed 
drainage flow at 1830, with stable conditions. This again 
Suggests that conditions on the south slope are different 


and isolated from conditions elsewhere. 
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4.3 Spectral Analysis of Winds 

During study of the wind speeds, there appeared to 
be a regular variation in speed, especially at the river 
valley slope stations. A simple correlation between station 
2 wind and the cross-valley wind component of station 3 
resulted in a correlation of 0.95, which suggested a 
periodicity since any random variation would be smoothed 
by the distance between stations. However, a periodicity 
does not seem probable as there is no apparent mechanism 
to result in such periodicities. The simple regression 
used is not suitable for the purpose of determining if the 
correlations reflect real periodicities as sensitivity is 
low. As an alternative the wind speed measurements were 


subjected to spectral analysis. 


A uly Ne ope CL bane : rogram 


The spectral analysis program used was devised by 
Lachapelle (1977) and uses the Fast Fourier Transform 
(FFT) as the basis for the modified periodogram method of 
analysis. The data were first preconditioned by the pro- 
gram as follows: 

1. The mean, variance and standard deviation were 

calculated. 

2. The mean was removed. 

3. Zeroes were added to the data set to reach a 
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4. The data were standardized by using the relation- 


ship 


Woe (X, OES 


where X is the mean, and S is the standard 

deviation. 
The standardized data were then passed to the FFT programs 
and the results made available in printed and computer- 
plotted form. Provision was made for smoothing the periodo- 
gram using a Daniell window. This is recommended by 
Kanasewich (1975), as it reduces the variance of the 
periodogram and produces a spectral estimate that will more 


likely be replicated if the experiment is repeated. 


Aes Cie cae Ging: 

Although the wind data were smoothed by the method 
of abstraction from the recorder charts, there was still 
a considerable amount of high frequency energy in the data 
which was of questionable value, as it added noise to the 
frequencies of interest. To remove these high Pree en cies! 
the data were filtered with a low-pass filter and as a 
result the characteristics of the filter used must be known 
to assess the effect on the data. The frequency response 
of the filter is the ratio of the amplitude of a wave of 
a given frequency in the series after filtering to the 


original amplitude before smoothing and this frequency 


aster en pad 


2 OGY Sarr g 33 wiht gt OWEted nat ors 


YES SU Re By 


# berated vi 


“Ob G44 EG aA gn haven Br) Loam BEM, ‘hava 


*at Wad 


“il behogmamos: PT, erat wai Vinthag a wag 
ond th wand fomy, ons auppiey af Pa i otanety ga twe 


ne it a“ 


tom J TEW Seas Bt erty 2a reat 18 aa * vooubog Bow nv re0be% 
aires | 1 RI 2 Aen 7 OARS . Wh": Aas - baausr asd pe 


. 
=) hs ay, i 


badiam, oid vd bedtoone aa hab: batw ada. fquaistA. 
cree cae SteRe cannes vate ne ond: sgn? aad 


ae fd i ee Vast \onmnipantt gre ta squons otdenabhs a 
i) 
ua? pd ge Pog botke 27 ae powh ee ofdenoFipaup: Yo. Bak : , | 


ietanewpay? Agi owens aveatgy at easat Wo ro hone pane 
$29 bas agate eeqc-wat' s Adie bevedht’. eee sh ad ; 


. 
Pe 


. 
edetd od Seam pany Hatt ty ont ie) eotzatveponreds ig He st 
astoges? vonneant ait “es as, i. toatts sid 126-04 
9p vas to abosttyma vite Oo o cae 


“at ot gatas, be 2 
, ‘ praeeet eva a ™ 


7 


| rica 
ere ft ail Poe \ i ot 


56 


response is a function of frequency.» The response function 


Offa fFiltes iwseEgieensby 
R(f) = f W(t) exp(-2rift)dt 


where f = + T = smoothing interval 


W(t) = weight function 


(4.1) 


This equation is the inverse Fourier transform of W(t) and 


therefore the inverse Fourier transform of the filtering 


function is the filter's response function. The absolute 


value of R(f) is 


[R(f)| = [(Re{R(#)})¢ + Giminr iat sc 
and the phase angle 


pentane Ink ray Renal 


The-aquiscrete. form sort the filter ais 


m 
Re) ean) We exp(i2afk) 
k=-n 


(4.4) 


Since it is desired that the filter not phase shift waves 
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at any frequency, Im{R(f)} must be required to be zero. 
This can be accomplished by requiring W(t) to be even 

(namely, W(-t) = W(t) or in the discrete case Cs Wh). 
If W(t) is even, the filter is called symmetric and the 


response function is 


R(f) = is Wit)imcos(emse)dt 


= OU le Wite a COS cai ty) dt (4.5) 
0 
or for discrete filtering having 2m + 1 weights 


R(f) Cos(ennk 
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Taking the simplest filter, an equally-weighted running 
mean or box car filter, where every weight is equal, the 


response function is computed from equation (4.5) using 
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Equation (4.5) gives 


RWries 2a oer ecostente) di 


SSS iat ) 


eR LE (4.7) 


The response function is shown in Fig. 4.10. The meaning 
of the negative response in Fig. 4.10 is that, as well as 
being changed in amplitude, the waves are reversed in 
poharyitys thateh,is, theaphasenisashifiteds180.degreess, qso 
that maxima are changed to minima and vice versa. This 
characteristic, plus the large value of the response varia- 
tion above the frequency of the first zero pole in the 
response (side-lobes), makes this type of filter undesirable. 


The Bartlett filter is defined as 


lcm l tlt A the | ties) 1/2 
W(t) = (4.8) 
0 Rule Te 


and the response function is 
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As can be seen (Fig. 4.11), the side-lobes are small, the 
first being 4.5 percent of unit response and all side-lobes 
are positive, indicating no phase reversal. This response 
is much more desirable in that it maintains phase relation- 
ships and is reasonably good at supressing frequencies above 
the first response pole. The Bartlett filter was chosen as 
the operating low-pass filter with the first zero pole set 
equal to the Nyquist frequency. 

In addition to the high frequencies, the trend must 
also be removed from the data, or the power contained in 
the trend will load the low frequency end of the spectrum 
and cascade or leak through to higher frequencies of the 
spectral estimate and invalidate the results of the power 
spectrum analysis at many frequencies. A high-pass filter 
was used for this purpose. The high-pass filter was synthe- 
sized from the low-pass Bartlett filter by first deriving 
the low frequency trend from the series, then subtracting 
this trend from the original time series. The resulting 
response function will be the inverse of that for the low- 


bass filter. 


(v8) 


sf) , Tere, oem cater sents ane” dd 
setol-ob la [ig baie we enone tne to 
cov ayy ,fstieyas ateng on pakt #3 


he 
roars ‘© 


“KOhIST 99 evade 2 estat he Ft J6ng ri gigsviadh jason 
ayoda setanecentt gnteewtque 34 baal ‘elannagney: a ar 
+s vowpas ee WRUTTT Goshinget ait aod aanogeay 
fae ofog aves Beek) wih tj cer he eshgodat® waraina 
» SOR steal ott on! ue Af 
ragh add J 2etoaorpart) eee ott ay aa? endian an 
sigkonsn young meh'oad poeahewng aaa’ novoniins 96 Dal : 
sioaqe vite te aso! y>ripenean) Met amg veal evi ind 8 
qi: to eotbabupa? swig hhad Hgiosnd seat vo sheseds 
ved att te eb taeen ott Shear abet me atnniian te 
coflh} gung-ughd A .cofcneepeet nee fe avrons 
-sitegd tow. yerhty\geea-iem oat “4 SETQTNG | ‘ahs i 
snivivsh S287 ya vast? rtehives. 2a8y-Wor eit 
Gal ue reat ‘yautae any an? sawn aha 
tie ro eaten mt Kaohg tro ete 
sun e-em Sa ‘a San ont: an rt te re 
| * fie a . Ps ee 


i 


a , 


is 


R( fF) 


60 


W(t) 


Vel 1/2 


Box—Car Filter 


Frequency 


Big® 4510¢%eBox-=car: fi Ltenvandmitsmresponse function: 


ery 1/2 
Bartlett Filter 


Frequency 


Fig. 4-Jjemebarclett tancer and 7ts response function. 


co 
Vis 


rf i¢ ri 


9 


ra 


tymy¥oqpandugsy eff twa oleae 


Hy 


6] 


4.333 “Confidence Limits 

To evaluate the time series, the significance level 
of the spectral peaks must be determined. Generally, 
confidence limits are generated by assuming a randomly- 
varying time series usually called white noise, that has 
a constant power per unit frequency. The test statistic 
is then the ratio of the sample variance of the line to the 
variance of the white noise. In many meteorological appli- 
cations, the use of a white noise source is not possible. 
One of the characteristics of many meteorological spectra 
is a general suppression of relative variance at high 
frequencies as compared to an even distribution of variance 
cChOSSminemspectrum expectedshyonea white noiseasource. 
This characteristic is caused by the appreciable auto- 
correlation found in meteorological data and has been 
called red noise (Gilman et al., 1963). The red noise 
spectrum may be derived by considering a first-order linear 
Markov process, where a discrete time series depends only 
on its own immediate past value plus a random component. 


This may be represented by 
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ot= pl) sethetpoputation! value of the 
lag 1 autocorrelation 


It has been shown by Gilman et al. (1963), that the 


estimate of spectral density as a ratio of red noise to 


white noise can be expressed in the form 


Dali) Se Balboa Bi 4.11 
5 5 (4.11) 
o-/m l= Zip COSw ON tip 


where uw) is the frequency 


eat is the variance of the red noise spectrum at 
frequency w 


scan is the variance of white noise per line 
m is the number of spectral lines 


In the spectral program, the spectrum is normalized, so 
that g? = | eenOneredanolse 10.) m= ail ie and if the sample 
has come from a parent population which represents the red 
noise process, then each sample autocorrelation value should 
provide the means to estimate p, the population lag 1 auto- 
correlation (Eddy, 1968). Since the values used should be 
those which are statistically significant and since the 


lag 3 autocorrelation value of the sample winds was around 


0.05, po was estimated from the sample by 
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where R(L) is the sample autocorrelation of lag L. 

The statistic to be tested during the red noise hypo- 
thesis is the ratio of the sample variance of the line to 
the red noise at that same frequency. For n pieces of 
independent data and m spectral lines, the equivalent number 
of degrees of freedom for the Chi-square distribution of 


a spectral line is 
Spey cole SN (4.13) 


Meteorological data are not independent due to the high 
autocorrelation, so that the number of degrees of freedom 


of the data is 
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spectral lines will follow the red Aeaee model. At those 
frequencies outside the pass-band of the filters, the 
spectral lines will contain far fewer degrees of freedom 
than those inside the pass-band, since the correlation of 
those frequencies outside the pass-band is much higher due 
to the filtering action. For simplicity, the assumption 
was made that all the degrees of freedom contained in the 


data would go into those spectral lines within the pass-band 
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where the red noise model was valid. The red noise model 
fit well in the spectrum of wind speeds from a frequency 
of 0.02 to 0.24 resulting in an effective number of 
spectral lines m' = 70 compared to 128 in the complete 


spectrum (Figs 74.12). 


Power 


Frequency 


Pigae4.\c. tbhrecte Or, pand-Dass Tiditering on spectrum 
compared to red noise model. 


The filtered spectrum, which is smoothed by a Daniell 
filter, has the number of independent spectral lines reduced 
by 1/Md, where Md is the number of points in the Daniell 
filter (for unsmoothed spectrum Md = 1). The degrees of 


freedom per spectral line is then 
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Confidence intervals were calculated using a Chi- 
square distribution. The confidence limits were obtained 
by multiplying the red noise value by x5 of and xe /v £OG 
the upper and lower limits, respectively, where a is the 
probability level selected. Since the degrees of freedom 
vary between one and two for the unsmoothed spectrum and 
since the distribution of the Chi-square function varies 
non-linearly between the values of v, graphs were plotted 
of the Chi-square function in order to obtain more accurate 
estimates of the confidence intervals for the unsmoothed 
spectrum. The same was done for the confidence intervals 
of the smoothed spectrum. 

The 95 percent and 5 percent confidence levels and 
the red noise spectrum were drawn on the graphs as shown 
in Fig. 4.13. All spectral lines exceeding the upper limit 
were tabulated and the significance levels were calculated. 
The confidence limits were also drawn on the filtered 
spectrum and the results compared to the unfiltered spec- 


trum (Fags 4 744,). 


4.3.4 Results 

Table 4.1 shows the results obtained from the uns 
smoothed periodogram and lists only the peaks for the prob- 
ability level > 95 percent and for percent variance 
explained greater than 1 percent, as peaks that explain 
such a small percentage of the variance would not be of 


much practical interest. Where there were two peaks to- 
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0.08 0.16 Un aeys 0.40 
FREQUENCY 


Unsmoothed periodogram analysis of station 5 
wind showing red noise spectrum and a 90% 
confidence interval. 


0.08 0.16 0.24 0.32 0.40 
FREQUENCY 


Three point Daniell-filtered periodogram 
analysis of station 5 wind showing red noise 
spectrum and a 90% confidence interval. 
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Periods (minutes) of significant peaks and 
significance levels (P %) for station winds 
from unfiltered periodogram. Cross-valley 
and along-valley wind components at station 
3 are indicated by V and U, respettively. 


P(%) 
% Variance 
99 Whe) 95 Explained 
Tale 0 4.6 
9.6 4.9 
Le3 4.6 
6.0 as 
[este 8.9 
Led Se 2 
Bale. 2 2AM 2 
19.2 725 
14.67 11.6 
19.2 8.0 


Indicates that the peak exceeded the 99.5% level. 
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gether near the same probability level, only one was listed 
as it is probable that adjacent peaks represent variance 

in the same frequency. With 70 spectral lines 3.5 would 

be expected to equal or exceed the 95 percent confidence 
level and, in all cases, this was observed, with 3 or 4 
peaks at or above 95 percent. Table 4.2 lists the signifi- 
cant peaks for the 3-point filtered spectrum. 

Inmelabier4s les tata onc: snows four significant peaks. 
although the peak at six minutes explains little variance. 
The peak at 7.9 minutes would seem to be the most signifi- 
cant and most likely to represent a real periodicity. How- 
ever a comparison with Table 4.2 shows that the peak is 
resolved at 8.1 minutes and has decreased in significance. 
This indicates that there is little likelihood of periodi- 
epties at’station 2. 

Neither wind component at station 3 shows a peak 
above 97.5 percent in Table 4.1 and in Table 4.2 little 
changes except that the peak at 18.1] minutes has moved up 
from a 95 percent to a 97.5 percent significance level. 
Thus it is doubtful that a real periodicity exists. 

Station o Shows two peaks in, japle 5.1, one at 14.6 
minutes explaining 11.6 percent of the variance and signifi- 
cant to better than 99.5 percent. This would appear to 
reflect a true periodicity, but chance could dictate this 
result.) INatathe periodicity ise real 1s) supported «in 
Table 4.2, where the period appears at 15.4 minutes and the 


significance level has increased to over 99.9 percent. 
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Table 4.2. Periods (minutes) of significant peaks and 
significance levels (P %) for station winds 
for 3-point Daniell-filtered periodogram. 
Cross-valley and along-valley wind components 
at station 3 are indicated by V and U, 
respectively. 


P(%) 

% Variance 

Station 99 We 95 Explained 
2 9.6 2°29 
8.1 20 
3-V Toon ba / 
(pele (bee 
3-U 19722 siefs 
5 15.47 | ee 


Indicates that the peak exceeded the 99.9% level. 
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The actual period was probably 15 minutes. Two questions 
must be answered. Does this periodicity make physical 
sense, and is there independent evidence of periodicity? 
The south slope is heavily wooded, which could act as a 
friction brake (reducing movement) until the air cooled 
enough to increase its density and thus its speed, where 
the increase in speed and the resulting mixing would slow 
down the cooling and thus slow down the air. As well 
station 5, unlike station 2, was located at the foot of the 
slope on a plateau, so that the measurements showed the 
effect of the entire slope on the drainage wind. There 
should have been a measurable effect on the temperatures 

at station 5, resulting in faster cooling during wind speed 
minima than when the wind speed waS a maximum. Unfortun- 
ately, the resolution of the thermograph was not sufficient 
to observe this effect and provide independent evidence 


of periodicity. 
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CHAPTER 5 


CARBON MONOXIDE 


5.1 Measurements 

Carbon monoxide (CO) concentration measurements were 
taken for a five minute duration at each monitoring station 
throughout the experiment (stations A-H; see Fig. 2.2). 

The resulting chart recording was divided into subsections 
and the CO concentrations abstracted from the chart by the 
Same method used for the wind speed chart recordings. The 
five minute series of measurements were then averaged and 

the resulting concentrations plotted in parts-per-million 

(ppm). Fig. 5.1 and 5.2 give the resulting CO concentra- 

tions from 1700 to 2400 MDT. 

Most concentrations, except for those stations close 
to River Road, were less than 1.0 ppm. Station A did not 
show any trend in concentration, but varied between 0.0 and 
0.8 ppm. Stations B to F showed a rising trend in CO 
concentration peaking around 2000 MDT. This was in opposi- 
tion to the pollutant release rate, as the average number 
of motor vehicles dropped from 1700 per hour at 1730 to 
900 per hour at 2000. This indicates that the increase in 
CO concentration due to the decrease in wind speed prior 
to the onset of the temperature inversion more than com- 


pensated for the emission rate decrease. The CO concentra- 
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tions at these stations dropped to average values that 

were lower after 2000 than before 2000. This was very 
evidenthatestetiongbuandtdramaticwat stationir whenetthe 
average concentration dropped from 2.1 ppm before 2000 to 
0.2 ppm after 2000. Station 6 showed the opposite trend 
and, although no measurements were taken at station H until 
after 2100, the concentrations were in accord with those 
measured at station 6 after 2100. These trends indicate 
the effect of the formation of the temperature inversion 
and the onset of drainage winds. 

The concentration at station C peaked at 0.8 ppm 
before 2000 and decreased by 2035 to 0.2 ppm. Station B 
showed a peak of 1.2 ppm at 1955, but the concentration was 
still 0.9 ppm at 2030, indicating that station B must have 
been within the inversion and came under the influence of 
the drainage wind later than experienced at station C. 
This is in agreement with the time taken for the inversion 
to deepen, as at 2145, when the inversion had filled the 
valley for one hour, the CO concentration at station B was 
near zero. The reason for the non-conformity of the time 
variation of CO concentration at station D has not been 
explained. Concentrations at station E showed the same 
Time Variationwas whose at stationic.. This is to.be 
expected, as both stations are close to the same height in 
the valley. The effect of the drainage flow on CO con- 
centrations was demonstrated most dramatically at stations 


F and G. These stations were on opposite sides of River 
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Road. The CO concentration variations were almost mirror 
images of one another as before 2000 the wind was southerly 
across River Road and after 2000 was northerly with the 
onset of the drainage winds. The rapid variation in con- 
centration at station G after 2100 was due to the close 
proximity of the station to River Road. At short distances 
single vehicles or groups of vehicles were registered on 
the CO meter. Station H was far enough away from the road 


for such variations to be smoothed out. 


s.a8 Potlucron Model 

As a quantitative check on the magnitude of the mea- 
sured CO concentrations, a simple Gaussian plume pollution 
model was used. For wind direction within 75 degrees to the 
perpendicular of a highway, Calder (1973) has shown that 
an approximate expression can be used to determine pollutant 


concentrations from an infinite line source. 
ee 
y neal or Gx) cos ¢ 


where X,, = the concentration at the perpendicular distance 
y y from the highway 


U = mean wind speed (m ar | 

Q = emission rate (g m7! eth) 

? = angle from the perpendicular to the highway 
) 


= the vertical standard deviation of plume 
distribution at the downwind distance x, from 
the source 
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The vertical standard deviation of plume distribution is 
dnaenpiricakgreliationsnips Choy Xe<e0sinkmrandsforsPRasquil | 
stability class F (moderately stable), Zimmerman (1975) 


gives the following expression for o (x): 


DJG) > Se AgOk & 05842)9-81558 


where X is in km. 


The emission rate Q can be found if the emission factor, 


] 


EF(g veh km!) and traffic volume, TV(veh bro!) are known. 


(EF) (TV) /3600 


yap) 
tl 


4 -| 


Dee TOdan CERN CVm(Geeat tkme 


As measurements of wind direction in relation to River Road 
were not available, a time was chosen when the wind speed 
measured at station 3 was a minimum. The assumption was 
made that, in the lower levels over the river valley floor, 
the wind speed equalled that of the drainage wind. This 
flow was perpendicular to the road so that ¢ = 0. The time 
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Environment lists the emission factor for the average mix 


Personal commumrcation wltneAr Schulz. P.End. 9 Of 
the Air Pollution Division. 
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of motor vehicles in Edmonton as EF = 23 g veh” | 


Traffic counts made by the City of Edmonton list the aver- 
age traffic volume at 2200 as TV = 500 veh hr7l, and this 
was confirmed by an independent count on the night of the 


experiment. Thus 


(2.78 x 107*)(23)(500) 


yan) 
iT} 


Shee Vena km! so! 


For comparison to the measurements which were in ppm, the 


calculated concentrations were converted to ppm by 


x 
x(ppm) = ane y 
a 


where o. = density of air 


4 


= Station pressure at YXD = 9.15 x 10) Pa 


T = temperature at 2200 = 275 k 
Ri segas constant — 267 J ka Ko! 
Because the vertical standard deviation of plume dis- 
tribution is a function of stability, concentrations were 
calculated using the equation for neutral stability to give 
an indication of the sensitivity of the model to stability. 


Zimmerman (1975) gives the following expression for o (x) for 


neutral Stabid ity: 
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a(x) = 34.459(X + 02722) ° 96974 


where X is in km. 
Using these figures, the results obtained are listed in 


Tra bsier:5s. 12, 


Table 5.1. Comparison between measured and calculated 
CO concentrations for moderately stable (MS) 
and neutral stability (N) conditions. 


Concentration (ppm) 
Calculated 


Station X(m) Measured MS N 
G 33 baad Ors Led! 
H 60 | ieee: 1.9 J igfsks 


The correspondence between the measured and calculated values 
indicates that the measured values were reasona le. An 
attempt to calculate CO concentrations before 2000 was not 
made because the angle of the wind to the perpendicular of 


River Road precluded the use of this simple model. 
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CHAPTER 6 


SUMMARY AND CONCLUSIONS 


During the evening of 21-22 October 1977, an experi- 
ment was run to obtain measurements of wind speed, tempera- 
ture and carbon monoxide concentrations in an effort to 
determine the micrometeorology of the North Saskatchewan 
River Valley as it passes through Edmonton, Alberta. The 
aim of this experiment and a series of others run and 
planned is to provide the information needed to create a 
simulation model of the microclimate and pollutant distri- 
bution within the river valley. In this experiment, 
number four of a series, six stations were instrumented 
On an approximate cross-section of the river valley near 
the center of Eamonton:. 

The results of the measurements showed that at 1800 
MDT the valley was isothermal up to the rim of the valley 
and slightly colder than the air above the plain of the 
city. This was thought to result from the earlier shading 
of the valley as compared to the plain. The prevailing 
wind as measured at the municipal airport and on the rim 
of the valley was from the southeast, and the wind over 
the river was easterly due to the channeling effect of the 
valley. The upslope-downslope wind component measured on 


the north (south facing) slope of the valley indicated up- 
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slope flow. The temperature at the foot of the south 
(north-facing) slope was 4°C lower than the rest of the 
valley due to the earlier onset of cooling of the south 
slope. The south slope drainage wind was well established 
by 1830, with no along-valley component evident. By 1930 
the temperature profile above the plain (urban profile) 
was isothermal and the lowest 25 m of the valley showed a 
slight inversion. The winds above the river were still 
easterly with a downslope wind component on the north slope. 
At 2000 an inversion existed in the lowest 14 m of the 
valley with a temperature gradient of 15°C/100 m. This 
deepened to fill the valley by 2050, with a gradient of 
5.2°C/100 m. The south slope was still colder than the 
rest of the valley and the resulting gradient, obtained by 
comparing the south slope station to the valley north rim 
station, was 13.9°C/100 m. An inversion had formed in the 
urban air by this time and was less intense than the valley 
inversion with a gradient of 3°C/100 m. During the time 
the inversion was filling the valley, the wind above the 
river had veered from easterly to become westerly, although 
the prevailing wind remained southeasterly. This was 
thought to be due to the geometry of the prevailing wind- 
river valley system, such that the westerly wind in the 
east-west section of the river valley was a compensating 
flow between two northerly flowing winds in north-south 
sections of the valley. After 2050 the wind above the rim 


of the valley backed to westerly from southeasterly, 
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apparently in association with the rising inversion. The 
Winds above the valley rim veered back to southeasterly 
by 2230 and the valley wind was still westerly. The down- 
Slope drainage flow was well established on both the north 
and south slopes of the valley, with the depth of the drain- 
age flow on the north slope at least three metres deep 
half-way up the slope. The urban inversion gradient was 
3.4°C/100 m and in the valley, as measured by the station 
on the river, the gradient was 6.0°C/100 m, while the 
gradient measured on the south side was 13.4°C/100 m. This 
situation remained stable until the conclusion of the 
experiment. At 2400 the winds were still westerly over 
the river and downslope on the slopes of the valley. The 
urban inversion was slightly stronger at 4.6°C/100 m. The 
inversion gradient over the river was 6.7°C/100 m while 
over the south slope the gradient was 15.1°C/100 m. 

Spectral analysis was performed on the slope winds 
and on the wind over the river. No evidence of periodicity 
was found for the valley wind or the north slope drainage 
wind. There was a periodicity of 15 minutes found in the 
south slope drainage wind. This was thought to be connected 
to the wooded nature of the south slope, but independent 
evidence for its existence was not found in the tempera- 
ture records. 

Carbon monoxide concentrations were measured at 
several points on the north side of the river (Fig. 2.2). 


The distribution was found to depend strongly on the wind 
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Field and was seen to decrease north of River Road and 
increase greatly to the south of River Road with the onset 
of the drainage wind. Over the flood plain of the river 
valley, the drainage wind was seen to extend to 15 m in 
depth by 2400 by observation of plumes emitted from the 
Glenora Club. A simple pollution dispersion model that 
assumed an infinite line source and a Gaussian distribution 
was applied to the vehicle emissions on River Road. Good 
agreement was found between the measured and theoretical 
Calculations for the case of airtlow perpendicular to River 
Road. 

The river station 3 was near the north side of the 
river and, in the last two hours of the experiment, it 
showed a northerly wind component. This, along with the 
drainage flow recorded on the slopes of the valley, implies 
that the air was converging over the river, and lends 
Support to the concept of a double-vortex circulation 
within the valley. It was not possible to look for the 
effect such a circulation would have on the temperature 
profile measured over the river due to the limitations of 
the equipment, which ceased functioning before the time 
any effect would be noticed. The lower temperatures on the 
south side of the valley could not be maintained with this 
type of circulation unless there was little interchange of 
air across the valley and a further study is required to 
resolve these points. 


The results of this study indicate that an inversion 
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stronger than the urban inversion can form in the valley 
and that such an inversion is accompanied by a valley wind 
system. This then has the potential to trap pollutants 
within the valley, especially on the south side where the 


inversion is strongest. 
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APPENDIX A 
TEMPERATURE DATA 
This appendix contains the temperatures for all 
stations used. This includes the Edmonton City Tower 


temperatures A, B, and C for 10 m, 45.7 m, and 91.4 m 


respectively. 
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Table Al. Temperatures at 10-minute intervals for valley 
stations used in vertical profile. 


Station 

Time ] 4 4a Ab 4c 4d 3 
1620 eZ es, 1032 howe ROre2 19023 

1630 T2380 VZ33 10.6 1025 10.4 TOmes 

1640 Te? breg 10.0 10.62 10.4 Toe) 

1950 Tie Ta's5 TOe5 102 L034 TORZ 

1700 Whee Va 10:25 TOs 3 105.2 VOas 

1710 ie Lia 10.0 1020 10.0 1070 

T2290 les Ue? 1022 9.9 9.9 9.9 10 
1230 1e0 a) LO 9.9 10.0 O70 10 
1740 10.8 1022 10.0 122 1032 HOw2Z 9 
1750 TORS less G4 10.0 9.8 9.9 9 
1800 9.9 1Oa93 9.6 Saar 9.6 9.6 9 
1810 10.0 10.0 G2 9.0 9.0 9.0 9 
1820 9.] 9.5 8.9 9.0 Bear g0 8 
1830 Sid 8.8 8.6 ad) 8.8 8.8 8 
1840 Sr 835 Stes) 8.6 S27 8.6 7 
1850 8.0 8x3 B22 oral 8.0 Ban i 
1900 Las Se 8.0 85.0 Vicks, lesh 7 
1910 Weak 8. ] ia Vals) 7.4 reds 7 
1920 7256 oe | 16 126 7.6 7.6 7 
1930 dish Wud Teal eee Wed roe) 6 
1940 Lee 136 Lave Eee dice Tae i 
1950 460 426 LO Timed yee) 6.9 6 
2000 6.9 7.4 ard 6.9 Ors (hed 5 
2010 640 7.3 636 6.6 6.10 6.6 4 
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Table A2. Temperatures at 10-minute intervals including 
Edmonton City Tower 


Station Tower Level 

Time 2 5 6 A B C 

1730 10.0 | 9.5 Bret 
1740 She) Os 8.7 
750 9.4 931 835 
1800 Ores 9150 8.4 
1810 8.9 orl 8.9 8.4 
1820 od! 8.6 8.6 8.1 
1830 8.4 S23 Sis 8.0 
1840 139 4.7 738 Oss: TET 
1850 Wad 4.4 Grace Sas 8.2 
1900 7.4 4.2 Be3 5 7.4 7.4 
1910 7.4 SR) Pee 1633 7.4 
1920 fel Sal 6.7 Peal Via ad 
1930 6.9 Sa Tt 6.4 ass Tee 
1940 Gia5 210 6's] 6597 Ts: 
1950 GPz3 a7) Sod 6.4 Te 
2000 9 220 G2 58 67,:3 Theek: 
2010 6.0 2.4 5a) OriZ 6n4/ 
2020 5130 244 5.6 On] G6 
2030 Diao 1.4 475 oe By. 6.8 
2040 4.3 ev] 57.0 one Pee 
2050 cE 0.8 Ge. Sra TiS 
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Tables A2. Continued: 

Station Tower Level 
Time s) 
2100 3.4 O55 4.4 4.4 5.0 
2010 Stl as: 4.4 Se0 
220 29 0.0 4.0 4.9 
2130 7a ooh -0.3 329 3.6 4.9 
2140 2D -0.6 3/26 4.5 
2150 2.4 -0.9 3.4 4. 
2200 Le -0.9 Sa 33 4.2 
2210 1.9 =e sec Shes 
2220 Tea7 -1.4 3761 328 
2230 1.4 -1.4 Bee Spal 39 
2240 Tel -1.4 Zo See 
2250 (had -1.6 2.4 36 
2300 ina! -1.9 ia Cie Sint 
2310 Tel -2.1 1.9 3.4 
2320 Om -2.2 agi sie 
2330 O73 -2.5 PAT Sil 
2340 O53 -2.8 1.4 37,0 
2350 0.0 -2.8 1.4 3.4 
2400 -0.3 -3.0 eee S00 
0010 -0.2 = Hes BGS. 
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